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The mechanism of delocalization of two-dimensional Dirac 
fermions with random mass is investigated, using a superfield 
representation. Although localization effects are very strong, 
one fermion component can delocalize due to the spontaneous 
breaking of a special supersymmetry of the model. The de- 
localized fermion has a non-singular density of states and is 
decribed by a diffusion propagator. Supersymmetry is re- 
stored if the mean of the random mass is sufficiently large. 
This is accompanied by a critical boson component. 

PACS numbers: 73.40.Hm, 71.55.Jv, 72.15.Rn, 73.20.Jc 



The transition from localized to delocalized states of 
non-interacting quantum particles in a random poten- 
tial is a phenomenon which is characterized by symme- 
tries. In contrast to classical critical phenomena, where 
symmetries are either discrete or compact continuous, it 
was observed that the transition from localized to de- 
localized states of a particle, described by a discrete 
Schrodinger equation (tight-binding model), are related 
to non-compact symmetry groups 0]. Nonlinear a models 
with the corresponding symmetries provide an effective 
large scale description, presenting the relevant degrees of 
freedom for localization or delocalization. They describe 
an effective diffusion of the quantum particle with diffu- 
sion coefficient D > 0. An important physical property 
of D in the nonlinear a model is its flow under renor- 
malization. In general there are fixed points, one for 
delocalized states (D > 0) and one for localized states 
(D = 0) In two-dimensional systems the renormal- 
ization always drives the diffusion coefficient to zero [Q, 
therefore reflecting the absence of delocalized states, at 
least in the absence of more complicated extensions of 
the Schrodinger equation like spin-orbit coupling. 

It turned out that for a number of interesting physical 
systems the effective quantum theory is not defined by 
Schrodinger particles but by Dirac fermions. The main 
reason for this is a linear dispersion and a sub-structure, 
either given by a sublattice or a spin. For instance, the 
integer quantum Hall transition (QHT) in a 2D electron 
gas with magnetic field can be formulated with Dirac 
fermions without a magnetic field |p|-[l0|| . Other examples 
for Dirac fermions are the degenerate semiconductor Jfi I 
and quasiparticles in a 2D d-wave superconductor |]l2| , |l3| . 

A Dirac fermion is a quantum particle with symmetry 
properties different from those of the Schrodinger parti- 
cles. In particular, the symmetry of the 2D Dirac Hamil- 
tonian is discrete in contrast to the continuous symme- 
tries of the Schrodinger Hamiltonian. This fact has im- 



portant consequences for the delocalization of the Dirac 
particle in d — 2 ||, and will be discussed in this letter. 
The Dirac Hamiltonian in 2D reads 



(1) 



Vj is the lattice difference operator in j-direction, M 
is the mass of the particle and Oj is a Pauli matrix. 
The localization properties of massless Dirac fermions 
with random vector potential was recently studied ||[l4j . 
It turned out that the low energy states are delocal- 
ized. Since the related Hamiltonian matrix has only 
off-diagonal elements this result can be compared with 
similar observations in ID systems: there are delocal- 
ized states at the band-center with a singular density of 
states (DOS) if the Hamiltonian represents hopping be- 
tween sublattices or different spin states |lq] . In contrast 
to this it of interest to consider models where also a di- 
agonal (potential) term appears in the Hamiltonian, and 
which have a non-singular DOS. 

The Dirac Hamiltonian Hd is an effective two-particle 
Hamiltonian because the Dirac theory includes particles 
and holes as the two components of the Dirac spinor. 
Hd is Hermitean and invariant under the transformation 
Hd — > —o?,Hdo~j,, provided the Dirac mass M is zero. 
However, this symmetry is not interesting here because 
it is always broken by the mass. Moreover, there is a 
space-dependent discrete transformation 



Hd — ► —SH'dS 



(2) 



for which the massive Hd is invariant. The 2x2 matrix 
S r is changing between o~i and 02 by going from one site 
to its nearest neighbor site, and H' D is obtained from Hd 
by a space- rotation of 7r/2 and a reflexion of the y-axis. 
(This is just an exchange of Vi and V2 in Hd-) 

In order to compare the Dirac Hamiltonian with the 
corresponding Hamiltonian H = V 2 + V of a Schrodinger 
particle in a random potential V , we extend the latter 
to Hs = (V 2 + V)<73. This Hamiltonian describes par- 
ticles and the corresponding holes, and can be used to 
express the two-particle Green's function for Anderson 
localization without a magnetic field. Hs is symmetric 
and invariant under a non-compact continuous symmetry 
under H s — > (ca + Si<ti + S2(J2)Hs{caa + S1C1 + s 2 o- 2 ) 
with the condition c 2 — s\ — s\ = 1. The role of the chem- 
ical potential in the case of Dirac particles is played by 
the Dirac mass, as it was earlier discussed by Ludwig et 
al. @. 

Transport properties can be evaluated from the two- 
particle Green's function B 
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K(r,r';e) = -(Tr 2 [G(r,r';ie)a 1 G 1 (r' , r; ie)^]), (3) 

where G(r, r'; ie) = (H + ie)~^> is the one-particle Green's 
function of Ho or H$, and (...) the average over random 
contributions in the Hamiltonian. For localized states the 
two-particle Green's function decays exponentially on the 
localization length. 

H$ is invariant under the transposition T of the matrix 
elements whereas the Dirac Hamiltonian Ho is not. It is 
convenient to write the two-particle Green's function as 
a functional integral 

G jr {r,r'-ie)G T klk {r',r^e) = 



Xr'j'Xrj^rk^r'k' exp(-S )X>*£>X- (4) 

Sq is a quadratic form of the four-component superfield 
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(5) 



with a complex component Xr and a Grassmann com- 
ponent \t r . The reason for introducing the superfield 
is that an extra normalization factor for the integral 
in Eq.([|) is avoided because of J exp(— St^V^Vx = 
det(iJ^ -(- ie)/ det(Ho + ie) = 1. It is crucial that Sq 
is not of the usual supersymmetric form Jl6[ | , where both 
diagonal elements are H + ie, if H T j£ H fT^ . This 
reflects a fundamental difference between the symmetric 
Schrodinger Hamiltonian H$ and the asymmetric Dirac 
Hamiltonian Ho for the construction of collective fields. 
In the following we will concentrate on the Dirac Hamil- 
tonian and refer to the literature for the case of the 
Schrodinger Hamiltonian [p] Jl^JT^ ]. 

In addition to the discrete symmetry of Ho there is an 
invariance of the action So for e — under supersymme- 
try transformation 
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with Grassmann variables ip and %j>. It is important to 
notice that the Dirac mass does not break this symme- 
try but only the term proportional to e in (^]), since U 2 
is not the unit matrix. Therefore, the field conjugate to 
the symmetry breaking field, which is quadratic in the 
superfield, must be studied in order to take the relevant 
symmetry properties into account. This leads to the col- 
lective field representation HkJ 



XrXr 



(7) 



The matrix elements Q r ,...,P r are 2x2 matrices, corre- 
sponding to the two components of ^S> r and Xr- 



Since the Dirac Hamiltonian Ho is usually obtained 
from a large scale (or low energy) approximation of a 
non-relativistic problem, there are several ways to intro- 
duce disorder which are motivated by the original con- 
densed matter systems. One starting point is, for in- 
stance, the network model of Chalker and Coddington 
p2] | for the QHT. This phenomenological description im- 
plies a random Dirac mass, a random energy term and 
a random vector potential ]l0| ]. Thus this model repre- 
sents a complex situation which also includes fluctuations 
of the magnetic field. Here we are interested only in the 
simplest possible case for the QHT of a system in a homo- 
geneous magnetic field. (Strong fluctuations of the ran- 
dom vector potential may drive the system into another 
universality class. At this point it is not clear if random- 
ness in the vector potential is relevant in the experiments 
on a 2D electron gas.) The QHT can also be described by 
a tight-binding model with a homogeneous magnetic field 
[^| in a random chemical potential. The latter would lead 
to a random Dirac mass. However, there was the argu- 
ment that the random Dirac mass alone does not present 
the generic situation for the QHT because the DOS is 
zero at low energy [^|, i.e. there are no bulk states even 
in the presence of disorder. It turned out though that 
these states exist if one goes beyond perturbation the- 
ory. This effect was also found in numerical calculations 
M,E0|. A consistent treatment of this nonperturbative 
contribution can be based on an effective field theory de- 
rived from the collective field Q |]] . This representation 
will be used in the following to discuss the breaking of 
the supersymmetry defined in (|^) and its consequences 
for the existence of delocalized states. 

Averaging over a Gaussian random Dirac mass M 
(where (M r ) — m and (M r M r i) — gS rr i) and transform- 
ing the functional integral to the collective field creates 
the new action J|,[l(| 



S' = - ^ Trg 4 (Q2) + lndetg(H + ie - 2rQr) (8) 

9 r 



with Ho = (Ho) and the 4x4 diagonal matrix r = 
((C3) 1 / 2 , (03) 1 / 2 ). Trg 4 and detg are the "supertrace' and 
the 'superdeterminant', respectively [ jl6| . In particular, 
the two-particle Green's function at r =/= r' then reads 

K(r, r'; e) = g- 2 {{@ r ,xi + ®r,2i){®v ,12 + 0^ 21 ))q. (9) 

The functional integral (■■■)q = J ... exp(— S')PQ can 
be approximated by a saddle point integration. A spe- 
cial saddle point is Qo = (m/4)7o — 1(77/2)73, where 
7j is the diagonal block matrix (<jj , o~j ) and 77 = irgp 
is proportional to the average DOS p ||. The sym- 
metry transformations are now applied to the saddle 
point solution tQot. The discrete transformation (|^) 
changes the sign of 77, i.e. the discrete symmetry of the 
massive Dirac Hamiltonian is spontaneously broken if 
lim £ ^o V 7^ 0. This is the case for — m c < m < m c 
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with m c = 2exp(— ir/g) 0. The supersymmetry trans- 
formation (^J), on the other hand, gives r\ — ► i]U 2 . Thus 
lim £ ^ V 7^ indicates also a spontaneously broken su- 
persymmetry. This behavior is analogous to a Heisenberg 
ferromagnet, where rj corresponds to the magnetization 
and e plays the role of the external magnetic field. How- 
ever, the situation is more complex for the Dirac fermions 
because the breaking of two symmetries is involved, a su- 
persymmetric and a discrete one. As a consequence of the 
supersymmetry of So at e — there is not just an isolated 
saddle point but a whole saddle point manifold, created 
by the symmetry transformation U. Therefore, the field 



Q' r = UrQoU^ 1 
= Qo iv 



T 7o 



— ilf> r ai —lpr'4'rO'3 



(10) 



controls the fluctuations around the saddle point man- 
ifold with U r = t*U t t and U~ 1 = tU t t* . U r is here 
the matrix U of Eq. @ in which the Grassmann vari- 
able ip is replaced by the Grassmann field i/i r . That 
means for the large scale properties the integration with 
respect to Q r can be restricted to an integration with 
respect to the field Q' r . Thus the critical (long-range) 
part of the random mass Dirac theory is controlled by 
the one-component fermion (Grassmann) field tp r . The 
bosonic (complex) field has only short-range correlations 
and, therefore, is localized by the disorder. The reason 
is that the bosonic field corresponds to the discrete sym- 
metry transformation (j^) which has a long-range mode 
only at the critical point where the order parameter rj 
vanishes. The latter is indeed the case because the local- 
ization length of Qn — Q22 and P%x — P22 increases like 
(m c — |m|) -1 / 2 as the critical value ±m c is approached 
from I to I < to c |J. This indicates a growing influence of 
these bosonic fields on the large scale properties. 

The expansion of (||) up to second order in the gradi- 
ents yields in general an action of the type fl|, |l^ , |l8| , [2lf 



: J d 2 rTrg 4 ( 73 Q' r ) + a J d 2 rTrg 4 (VQ' r • VQ' r 



d 2 r 



5^ e/w T>g 4 (Q , r V M Q' f .V 1/ Q / r ), 

//./,' 



(11) 



where e M „ is the antisymmetric unit tensor, and the pa- 
rameters a and /3 are determined by the model. In par- 
ticular, for the quantum Hall effect there is a = <j X x, the 
(unrenormalized) longitudinal conductivity, and (3 = o~ xy , 
the (unrenormalized) Hall conductivity |2l|| . The topo- 
logical term / cPr£^e /1I ,Tirg 4 (Q / r V /t Q / r V,,Q / ,.) takes 
care of the Hall plateaux because the latter are a conse- 
quence of the (topological) edges states in the presence 
of localized bulk states. At the QHT, however, transport 
is dominated by delocalized bulk states. Therefore, the 
topological term should not play a crucial role in this 
case. In fact, for the Dirac Hamiltonian Hp with to = 0, 
i.e. for the choice Q' of Eq.rtnT), the topological term 



vanishes. The only terms which remain in the action are 
the linear off-diagonal elements of Q' 



S"={l/-wp) / d 2 rVv(e + £>V 2 )Vv, 



(12) 



where the average DOS p and the diffusion coefficent D 
can be evaluated form the saddle point equation. This 
surprisingly simple result, which satisfies the Ward iden- 
tity K(q = 0, e) = irp/e, reflects the fact that only a one- 
component Grassmann field contributes to the massless 
fluctuations, created by the broken supersymmetry. That 
means there is a simple physical structure for the well- 
delocalized Dirac fermions in the vicinity of to = 0. The 
divergent localization length of two real boson compo- 
nents will eventually turn into a restoration of the super- 
symmetry, where lim e _^o fj — 0. Since the supersymmet- 
ric theory in 2D does not have delocalized states |l6| , the 
restoration of the supersymmetry must be accompanied 
by a transition into a localized regime. This is the regime 
characterized by the Hall plateaux. The critical behavior 
of the two real fields (which can be considered as the two 
components of one complex boson field) at m = ±m c 
due to the spontaneously broken discrete symmetry (||) 
invalidates S" near these points. It must be replaced by 
a more complicated field theory which includes both the 
critical boson field and the critical Grassmann field of 
(O). This would require an additional matrix field 



Qr@~ 3 







-ip r a 3 



(13) 



added to Q' r in (Q . The real field components q r and p r 
are related to Qii, r — Qii,r and Pu, r — P22,r, respectively. 

As a direct consequence of these results the value of 
the conductivity at to = (the "conduction peak" ) can 
be evaluated from the Einstein relation a xx = (e 2 /K)Dp 



1 



■nh 1 + g/2-K' 



(14) 



This is in agreement with experimental results |2jj2^] and 
other theoretical work [§6],^7| ■ F° r weak disorder the sec- 
ond factor can be neglected. In this case the peak value 



is just the universal constant a x 



2 /irh. The latter 



was obtained for Dirac fermions in a random vector po- 
tential |5) and for the lowest Landau level with random 
spin scattering [28), regardless of the strength g. Thus 
the extra factor in Eq.(|l4]) indicates that the random 
Dirac mass, representing a random scalar potential, has 
a stronger effect on the Dirac particle than the random 
vector potential or a random spin scattering. 

In conclusion, the large scale behavior of the two- 
particle Green's function of 2D Dirac fermions with ran- 
dom mass is characterized by a single massless Grass- 
mann field on the interval [— to c , m c ] of the average Dirac 
to. It describes delocalized states with non-singular DOS 
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due to a broken supersymmetry by a diffusion propaga- 415 (1993) 

tor. There is a two-component real bosonic field which 
has a divergent localization length as the critical points 
±m c are approached. It corresponds to a broken discrete 
symmetry of the Dirac Hamiltonian. This mechanism of 
derealization is different from the one which is respons- 
able for delocalized states in the random vector potential 
or for random spin scattering. 
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